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Abstract. - We introduce an efficient method to solve the Bose-Hubbard model. The 
Schrodinger equation is solved by the successive construction of doorway states. The ground 
state wavefunction derived by this method contains all relevant many-body correlations intro- 
duced by the hamiltonian, but the dimensionality of the Hilbert space is greatly reduced. We 
apply the doorway method to obtain the chemical potential, the on-site ffuctuations and the 
visibility of the interference pattern arising from atoms in a one-dimensional periodic lattice. 
Excellent qualitative agreement with exact numerical calculations as well as recent experimental 
observations is found. 



The physics of strongly interacting quantum systems has been subject of intense investi- 
gations since the early days of quantum mechanics. With the recent developments in the field 
of ultracold quantum gases, systems of strongly correlated bosonic and fermionic atoms with 
tunable interactions have become amendable to experimental studies [1-5]. Complete theo- 
retical understanding of these coupled many-body systems is difficult, and only rare examples 
exist, for which analytical solutions for the ground state could be explicitly been given [6]. 
The description generally has to rely upon approximation methods such as the density-group- 
renormalization [7], a Gutzwiller ansatz [8-10] and quantum Monte Carlo approaches [11]. 
Of particular importance is the Bose-Hubbard model, which serves as a prototype system 
exhibiting a quantum phase transition [12]. Triggered by the recent observation of this phase 
transition from a superfluid to a Mott-insulating phase with bosonic atoms in one- (ID) 
and three-dimensional (3D) optical lattices [3,4,13], the ground state and correlation prop- 
erties of the Bose-Hubbard model have been extensively investigated by a large number of 
groups [14-16]. 

In this Letter, we explore an alternative and particularly efficient theoretical approach 
to solve the Bose-Hubbard model. The method is based on the construction of doorway 
projectors as initially introduced by Feshbach in the context of nuclear scattering [17]. A 
similar approach appears in some other areas of physics under different names. For example, in 
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solid-state physics it is known as the Lanczos method [18], in quantum optics as the resolvent 
method [19] and in physical chemistry as the Feshbach-Fano partitioning method [20]. In 
principle, any state contained in the Hilbert space of the Hamiltonian can serve as a starting 
state for the doorway method. However, the efficiency of the method can be greatly increased 
by choosing a start state which is contained in an appropriate subspace determined by the 
intrinsic symmetries of the hamiltonian. The doorway method thus genuinely leads towards 
the full many-body wavefunction of the system and it significantly reduces the calculational 
effort, yet maintaining all correlations of the many-body states. 

In this paper, we use the symmetry of the ground-state of the Bose-Hubbard hamiltonian 
to choose the starting state and to generate the correspondent sequence of doorway states 
leading into the full wavefunction of the ground state. In order to demonstrate the applicability 
and the advantages of this method, we determine the energy spectrum of the ground state, 
which features the phase transition, in both the superfluid and Mott-insulator regime in the 
homogeneous case. Furthermore, the correlation properties of a bosonic gas in an optical 
lattice are calculated over the full range of coupling parameters and qualitatively reproduces 
the measurements of [13,21], even though these experiments are performed in inhomogeneous 
systems due to the presence of an external trap. 

The doorway method consists in splitting the full Hilbert space of the many-body hamil- 
tonian H into different subspaces which are accessed by the action of projectors on H. Let 
P describe a projector onto a subspace P. For the sake of simplicity, we assume the initial 
subspace to be one-dimensional with P = \p){p\. Note, however, that P does not necessarily 
need to be one-dimensional. The first doorway state is given by \d^^^) = — P)H\p) 

where M^^^ denotes the normalization, the second doorway by the normalized basis state 
jd*^^)) = jp^i'^ — P ~ Di)H\S-^^) and so forth. The procedure terminates when (1 — P — Di — 
...-L>„)if|d("-^)) = with Di = (i = 1, 2, . . . , n) denoting the projectors onto the 

doorway spaces. The doorway Lanczos, resolvent or Feshbach-Fano partitioning methods have 
the same algorithm to generate the so-called doorway states. However, the crucial point is the 
choice of the most adequate starting state \p) to generate in an efficient way the solution. The 
space P® Di ® • • • ®P„ then represents an irreducible subspace for the Hamiltonian containing 
the initial space P. The many-body Schrodinger equation is solved in this subspace 

[E - Hpp - Hpn, p „ ^ —Hn,p)Pm = 0, (1) 

where Hpp = PHP, HoiDj = DiHDj, HoiP = DiHP with i,j= 1, 2 ... n. Recursively, one 
thus obtains the diagonalization of the Hamiltonian in this subspace. All orders of interactions 
described by the Hamiltonian are taken into account inside this subspace. Therefore, the 
expansion is not perturbative in the couplings, but it represents the full solution in a subspace 
of lower (reduced) dimension. 

We apply the doorway method to a system of N bosons in a ID optical lattice of M 
sites in the homogeneous case. The system is described by the Bose-Hubbard model H = 
^J'^Z'jLii'^'j^j+i + h.c.) + Y Ylj ^jin-j ^1)- which can be seen as a many-mode approximation 
for an ultracold gas of bosonic atoms occupying the lowest band in a periodic array of potential 
wells [9.12]. The operators a] and Uj are the creation and annihilation operators, respectively, 
of a boson in the j-th site, and nj is the number of atoms at the lattice site j {j = 1 . . . M). 
The strength of the on site repulsion of two atoms on the lattice site i is characterized by U. 
The parameter J is the hopping matrix element between adjacent sites which characterizes the 
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strength the tunnehng term. A periodic boundary condition is included aUowing for tunneling 
between the first and the last site of the lattice (site M + 1 = site 1). 

The Bose-Hubbard hamiltonian with periodic boundary conditions conserves two impor- 
tant quantities: the total number of particles and the total quasi- momentum. As a conse- 
quence, its Hilbert space can be partitioned into M decoupled subspaces corresponding to the 
values of quasi-momentum (q=0,...,M-l). Furthermore, the subspace with q=0 and, for even 
M, also the subspace with quasi-momentum q=M/2, can be further partitioned into symmet- 
ric and antisymmetric parts (" odd-even" symmetry) [22] . The ground state of the system lies 
in the symmetric part of the q = subspace. Thus, the initial P subspace is chosen within this 
manifold. The full set of doorway states constructed from this initial state will finally span 
the symmetric part of the subspace with zero quasi-momentum, providing the full solution for 
the ground state. In addition, the ground state of the system has two different limiting cases: 
the superfluid phase for J ^ U where the particles are delocalized over the whole lattice, and 
the Mott-insulating phase for J ^ U where the particles are localized with an equal number 
of atoms per site uq = N/M. 

For the commensurable situation (N = tiqM) we choose \p) ~ Ino)***^ corresponding to 
the ground state of the system in the pure Mott case {U ^ J). If an extra particle is added to 
or removed from the system {N ± 1 = uqM ± 1). The appropriate start state is chosen from 
the symmetric subspace, i.e., \p) = -;^J2jLi K'^o i l)jnf'^). Starting from these states, 
the doorway states are constructed following the method described above. Depending on 
the parameters of interest, one may also start from the pure superfluid state (J U) as 
the initial state in the construction of the doorways and follow the same procedure as before. 
Thus, the energy spectrum and eigenvectors are obtained by solving Eq. as a function of the 

parameters J and U . One can check that the superfluid ground state \'iPsf) ^ (X^fii '^]) I*-*) is 
asymptotically reached for J/U — > oo. The doorway method genuinely reduces the numerical 
effort to the diagonalization of the very small, but relevant part of the total Hilbert space. 
A good measure for the convergence of the method is given by the "reliability" defined 

as 7^ = \{iIjo\4'o°°^"'°'^)\ — X]r=o Q^iKV'oMi)!^, where |V'o) being the exact ground state of the 
system obtained from the exact numerical diagonalization of the Bose-Hubbard Hamiltonian, 
^^doorways^ being the ground state of the system obtained from the doorway method and n 
being the maximum number of doorways to generate the truncated solution. The meaning 
of the reliability can be understood in terms of the portion of the full symmetric subspace 
which has been covered by the solution. The convergency of the method scales linearly with 
the number of particles for the commensurate and incommensurate cases. Fig^a) illustrates 
this convergency for M = N — 5, considering the truncated solution for n = 2,3,4 and 5 
doorways. Excellent convergency is already achieved after the fifth doorway jds) is created, 
showing the subspace spanned by the doorways already include the most relevant portion 
of the full symmetric subspace. For this particular case, the Hilbert space is of dimension 
{^~^N~^) ~ ^'^^ whereas the symmetric zero quasi-momentum subspace has only dimension 
16 and includes the exact ground-state solution. The convergency of the method is very 
fast, since the recursive construction of the doorway states "automatically" creates the most 
efficient path into the irreducible subspace depending on the specific values of the parameters 
U and J. In fact, by applying the doorway method for larger numbers of particles and sites, we 
find that the number of doorways needed to achieve convergence roughly equals the number 
of particles. 

We derive the ground state energy and the chemical potential (pn/U ~ [E^^i — Ejq]/U) 
as a function of the dimensionless tunneling parameter J/U. Fig. Qlb) shows the chemical 
potential as a function of J/U comparing the exact solution of the Bose-Hubbard hamiltonian 
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Fig. 1 - (a) Reliability as a function of the tunneling parameter for N=M=5 and the maximum 
number of doorways n = 2,3,4,5. (b)Chemical potential fiN /U as a function of the dimensionless 
tunneling parameter J/U for M = 5 lattice sites and A'^ particles. The solid lines are the result of the 
doorway method using seven doorways. The dashed curves represent an exact numerical solution of 
the Bose-Hubbard model. 



(dashed lines) with the result of the doorway method (solid lines) for the case of M = 5 sites. 
The distance between the apparent parallel lines of Fig.^b) is almost 1/M, which vanishes 
in the thermodynamic limit (A'', M ^ co and N/M finite), characterizing the phase transition 
point where the two limiting curves ^uqAI and UngM+i are the same. For convenience, the 
doorway method has only been applied to the limiting cases N — uqM and N = tiqM ± 1 
where the determination of the corresponding initial states is described above. The curves 
inside the lobes are only shown for the exact calculation and they represent the evolution of 
the M-degenerate states of J = as a function of the hopping parameter. 

As a second application of the doorway method, we determine the on-site number fluctu- 
ation of the ground state cr,; = J (nf) — {riif' as a function of U / J and as a function of J/U. 
The results are shown in Fig. |2| Depending on the regime, the most efficient construction 
of the doorways starts from either a pure superfluid state (left graph in Fig. [21 or the Mott 
state (right graph in Fig. |2J|. For small C//J, the on-site fluctuations are found to decrease 

linearly starting from the initial value given by ^-^-(1 — -p). At intermediate U / J, the de- 
crease curve changes slope in the region of the phase transition and finally approaches the 
dependence 2^y^no{no + 1), which is identical to the perturbative result for the ground state 

wavefunction | V') — \no)^^^ ~^TjJ2{i j) ol'^i I'^o)'**^ with denoting the next neighbor sites. 
The results of the doorway method agree well with an analysis by Roth and Burnett based 
on the calculation of the superfluid fraction under twisted boundary conditions [14]. Starting 
from the pure Mott-insulator state, we derive the on-site number fluctuations for J/U <^ 1 
(right graph in Fig. |5J. In the pure Mott phase, no correlations among the particles exits 
and the one-site number fluctuations vanish. As J/U increases, correlations are created re- 
sulting in a steady increase of the fluctuations. After the phase transition region is crossed, 
the fluctuations approach the asymptotic value for U/ J = 0. It is instructive to compare the 
result with the commonly used Gutzwiller ansatz [9, 10] which, by construction, does not 
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Fig. 2 - On-site number fluctuations calculated for Al — 6 sites and TV = 6 particles using eight 
doorways to construct the ground state of the system. The initial state in the recursion is either the 
pure superfluid state (left graph) or the pure Mott-insulator state (right graph). The dashed line 
gives the result of perturbation theory for J/U <^ 1. 



provide particle correlations and thus genuinely underestimates on-site fluctuations at small 
J/U. The short-range correlations of the particles have to be included by additional pertuba- 
tive terms [10, 16]. Since the groundstate wavefunction constructed by the doorway method 
contains all relevant many-body correlations, the on-site fluctuations are rendered correctly 
for small J/U as well as for the region of the phase transition at intermediate values of the 
tunneling parameter. 

To demonstrate the efficiency of the doorway method for the determination of important 
experimental observables, the interference of atoms released from an optical lattice is modeled 
and the dependence of the visibility and the width of the first resonance peak of the matter- 
wave interference pattern on the parameter J/U is derived. Comprehensive studies of these 
two quantities have recently been performed by Gerbier et al. [21,24] for 3D optical lattice and 
by Stoferle et al. as well as M. Kohl et al. for ID optical lattice [13]. The observed interference 
pattern is directly related to the Fourier transform of the single-particle density matrix pij = 
(ajoi)- Therefore, measurements of the visibility are particularly sensitive to short-range 
coherence induced by the tunneling of the particles. Following the derivation in Refs. [14,25], 
the intensity is calculated for a one-dimensional lattice as I = Sfj=i "^''^P ~ '^j)) 

where k denotes the wave vector of the expanding matter wave after release from the lattice 
point Ti or rj , respectively. The result of the doorway method for N = 8 particles distributed 
over M = 8 sites, starting from the pure Mott state and taking nine doorways into account, is 
shown in Fig. |3 Periodic boundary conditions assure that limited size effects are minimized. 
Fig. Ota) depicts the resulting intensity distribution for different tunneling parameters. One 
clearly sees the development of sharp multiple-beam interference maxima with increasing J/U, 
corresponding to a decrease of the well depth in the experiments. For very small J/U, the 
remaining weak modulation of the intensity distribution indicates small-scale correlations even 
deep in the Mott-insulator regime, as have been observed experimentally [13,21,24]. The 
intensity pattern for small J/U is well fitted by the result of perturbation theory for an infinite 
number of sites given by / = 1 + 8^ cos(fcd) (dashed line in Fig. 13) where d is the separation 
of two adjacent lattice sites. 
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In Fig.|3|^b) the visibility, defined slsV — (/max — ^min)/(-^max + ^min), is plotted over a wide 
range of the parameter U/ J. The power of the doorway method becomes apparent by the fact 
that all features of the visibility are correctly reproduced over the full range of U/ J by using 
a ground state wavefunction created by the superposition of few doorway states including the 
initial Mott-insulating state. For large U/J (Mott-insulator regime), the visibility approaches 
the non-vanishing perturbative result V = 4(no + l)J/U [21,24,26] reflecting the persistence 
of short-range correlations. In the ballpark of the phase transition the visibility changes its 
analytical dependence on J/U and finally reaches a value close to unity for small U/J. This 
graph qualitatively describes the experimental findings of Gerbier et al. [21], even though the 
present model neglects the external trap and it is not ID. A detailed investigation of the 
visibility in a 3D- inhomogeneous case can be found at ref. [27]. In addition, we can directly 
calculate the width of the first interference peak for different number of sites and compare it 
with the experimental measurements of refs. [13] for ID optical lattice. In Fig.lSJc), the width 
of the first interference peak is calculated using the doorway method for different number of 
sites (M = 3, 4, 5, 6, 7, 8). We observe that the peak width starts to increase at a critical value 
close to U/J ^ 4 which characterizes the phase transition region. Our calculation describes 
qualitatively the results of [13] close to the phase transition point . At larger U/J > 10 the 
width of the interference peak in our model is limited by the finite number of sites. 

In conclusion, we introduce an alternative and powerful theoretical method to solve the 
Bose-Hubbard model for interacting bosonic particle in a periodic lattice. The method is based 
on the successive construction of doorway states which genuinely takes profit of symmetries 
included in the hamiltonian. It naturally terminates when the full irreducible subspace con- 
taining the starting state is spanned by the doorway states. The intermediate doorway states 
depend on the interaction and tunneling parameters of Bose-Hubbard hamiltonian providing 
the closest approximation to the ground state. Therefore, the convergence of the method, 
as can be quantified by the "reliability" introduced in this paper, scales as the number of 
particles. All relevant correlations of the many-body system are included in the ground state 
wavefunction derived in this way, as was shown by calculating the energy spectrum, the on-site 
fluctuations and the expectation value for single-particle correlations. We have limited our 
discussion on a one-dimensional model to keep the calculational effort on a low level. Concep- 
tually, there are no limitations for extending the method to 2D, 3D and to the inhomogeneous 
case. For a more detailed comparison with experimental flndings, in particular concerning 
the effects of an additional trapping potential [27], we will also add a harmonic term in the 
Hamiltonian leading to the appearance of different phases within the sample. Other impor- 
tant observables currently under experimental investigation, such as two-particle correlations 
of Hanbury-Brown and Twiss type [28] , can also be approached by the doorway method. Since 
the doorway method allows direct access to the many-body wavefunction of the system, it can 
be applied to other important systems involving many-body correlations such as fermionic 
gases coupled by Feshbach resonances or the fermionic variant of the Bose-Hubbard model. 
In addition, dynamical properties of the Bose-Hubbard model, such as many-body tunneling 
rates, may be investigated. 
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Fig. 3 - (a) Interference pattern of a matter waves released from a one- dimensional lattice for different 
values of the tunneling parameter J/U for eight particles distributed over eight lattice sites. The 
intensity is plotted as a function of the accumulated difference in phase after the expansion, (b) 
Visibility of the interference pattern as a function of the inverse tunneling parameter U / J. (c) The 
width of the first interference peak as a function of U /J for different number of sites M = 4, 5, 6, 7, 8. 
The dashed lines in (a) and (b) show the result of perturbation theory for an infinite number of sites 
but finite average population per site no = 1. 
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